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In this paper we generalize Yu’s [Ann. Probab. 24 (1996) 2079-2097] 
strong invariance principle for associated sequences to the multi¬ 
parameter case, under the assumption that the covariance coefficient 
u{n) decays exponentially as n ^ oo. The main tools that we use 
are the following: the Berkes and Morrow [Z. Wahrsch. Venn. Gebi- 
ete 57 (1981) 15-37] multi-parameter blocking technique, the Csorgd 
and Revesz [Z. Wahrsch. Verw. Gebiete 31 (1975) 255-260] quan¬ 
tile transform method and the Bulinskl [Theory Probab. Appl. 40 
(1995) 136-144] rate of convergence in the CLT. 


1. Introduction. Among various concepts introduced to measure the de¬ 
pendence between random variables, association deserves a special place 
because of its numerous applications and its relatively easy mathematical 
manipulation. A finite collection (Xi,... ,Xm) of random variables is said 
to be associated (or satisfies the FKG inequalities) if for any coordinatewise 
nondecreasing functions /, g on R™', cov(/(Ai,..., , A,„)) > 0, 

whenever the covariance is dehned. An infinite collection of random variables 
is associated if every finite sub-collection is associated. This concept was for¬ 
mally introduced by Esary, Proschan and Walkup (1967), who also deduced 
some of its most important properties. 

In the past few decades, a lot of effort has been dedicated to prove limit 
theorems for random fields {Xj)j^ 2 ,d- of associated random variables. In the 
case d=l, this culminated with the strong invariance principle of Yu (1996), 
from which one can easily deduce all the other major limit theorems, like the 
weak invariance principle and the functional law of the iterated logarithm 
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(FLIL). The present paper was motivated by the need for a similar result 
in the case d>2, which arises in the context of higher-dimensional models, 
like the percolation model of Cox and Grimmett (1984). 

The first asymptotic result for zero-mean associated random fields was the 
central limit theorem (CLT) proved by Newman (1980) for the (strongly) 
stationary case. This result says that if the finite susceptibility assumption 
holds, that is, cr^ := P(0 < OO; where p{j — k) := cov{Xj,Xk), then 

( 1 ) n-'^/^Sn^N{0,a^), 

where Sn ■= J2ji<n '' 'J2j^<n^j- This was generalized by Cox and Grimmett 
(1984) to the nonstationary case, under the assumption that u{n) —> 0 as 
n ^ oo, where 

(2) u(n) := sup ^ cov{Xj,Xk) 

k : IIj — 

and ||i|| := maxs=i^,,,^rf |zs|. 

The weak invariance principle for (strongly) stationary associated random 
fields was proved by Newman and Wright (1981, 1982) in the case d = l and 
d = 2, under the same finite susceptibility assumption. They also conjectured 
that the same principle holds for d> 2. A partial solution to this problem was 
given by Burton and Kim (1988) in the stationary case, and by Kim (1996) 
in the nonstationary case, under the finite r-susceptibility assumption: 

where Sn := T,j<N and [A^] := nf=i N's for iV = (iVi,..., Nfi) £ Z^. (If i,j G 
Z'^, we use the notation i < j if is < js, Vs = 1,..., d and i < j if is < 
js, Vs = l,...,d.) 

The conjecture was fully solved by Bulinski and Keane (1996), who proved 
that for a zero-mean (weakly) stationary associated random field iXj)j^ 2 ,<^ 
with uniformly bounded moments of order s > 2 and a power decay rate for 
the covariance coefficient u{n), we have 

(3) W4-)^W{-) 

where Wn{t) ■■= ^ = (W^(0)te[o,i]‘* is a d- 

parameter Wiener process with variance . We note, in passing, that for 
d = 1, generalizations to the nonstationary case and to the vector-valued 
case were given by Birkel (1988a) and by Burton, Dabrowski and Dehling 
(1988), respectively. 

The FLIL for associated sequences was obtained by Dabrowski (1985), 
under the finite r-susceptibility assumption with r = 1 and a condition which 
requires that E(Slfi)ln converges to with a power decay rate. 
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The strong invariance principle proved by Yu (1996) strengthened and 
unified all of these results in the case d = \ and implied other asymptotic 
fluctuation results, like the Chung’s type of FLIL for the maxima of par¬ 
tial sums; see Theorems A-E of Philipp and Stout (1975). More precisely, 
Yu showed that if [Xj)j^z+ is a sequence of associated random variables 
such that the moments of order s > 2 are uniformly bounded, the vari¬ 
ances are bounded below away from 0 and the covariance coefficient u{n) 
decays exponentially as n —> oo, then it is possible to redefine the origi¬ 
nal sequence on a richer probability space together with a standard Wiener 
process W = (lY(t))ig[o^oo) such that, for some e > 0 

Sn-W {al) = 0 (n^/^"^) a.s., 

where cr^ := E{S‘^). As far as we know, there are no generalizations of this 
principle to the case d>2. The purpose of the present paper is to fill this 
gap and to provide a powerful approximation tool that can be used in higher 
dimensions. 

Unlike the case d=l, the strong invariance principle for associated ran¬ 
dom fields in higher dimensions holds only for points N G Z'j[_ which are not 
“too close” to the coordinate axes. This is not at all surprising and a similar 
fact happens for mixing random fields; see Theorem 1 of Berkes and Morrow 
(1981). The reason for this phenomenon is the irregular behavior of E{S‘^) 
close to the coordinate planes. 

We proceed now to introduce the notation that will be used throughout 
this paper. 

Let (Yj)jg 2 '( be a weakly stationary associated random field with zero 

mean and p{j — k) := E{XjXk), y j,k G Z'^. Let u{n) be the covariance coef¬ 
ficient defined by (2). Because of stationarity, we have u{n) = Z]iGZ<^:||i||>n 
for every n > 0. We will suppose that / 9 ( 0 ) > 0 and cr^ ;= u(0) = P(*) < 

oo. 

For any finite subset V C Z^|_, we let \V\ be the cardinality of V, S{V) ;= 
Ejw Yj,cj2(U) :=E[S'^{V)] and Eyix) := P{S{V)/aiV) < x),x eTi. Note 
that for any finite subset V C Z^|_, 


( 4 ) 


K0)< 


14^1 


<cr^ 


Most of the time we will work with “rectangles” V C Z'^ of the form V := 
(a, 6 ] = with as,bs G Z 4 . U {0},as < 6 ^; note that |U| = [b — a]. 

We denote with A the class of all the subsets V of this form. 

We will use the following conditions: 

(Cl) supjg 2 ‘' < 00 for some r, <5 > 0. 

(C 2 ) u{n) = 0(e“'''") for some A > 0. 
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(C2') u{n) = 0{n for some u >0. 

We recall that a d-parameter Wiener proeess W = {Wt]t G [0,oo)'^} with 
variance is a Gaussian process with independent increments such that 
W{R) has a A^(0,cj^|i?|)-distribution for any rectangle R (|i?| denotes the 
volume of R). Using the same notation as Berkes and Morrow (1981), we 
put 


G'r Pi £ Z([.: js > rG(0,l). 

s=l I s'^s ) 

Here is the main result of this paper. 

Theorem 1.1. Let d > 2,r G (0,1) and (Wj) be a weakly stationary 
associated random field with zero mean and p{j — k) := E{XjXk) for any 
j, k G Z^. Suppose that p{0) > 0 and := pif) < oo- 

If (Cl) and (C2) hold, then without changing its distribution, we ean 
redefine the random field ® richer probability space together 

with a d-parameter Wiener process {Wt,t^ [0,oo)'^} with variance such 
that 


Sn -Wn = 0{[N]^/‘^-^) a.s. 

for N G Gr- Here e is a positive eonstant depending on the field {Xfi-^r^d . 

From the previous theorem one can easily deduce the following CLT: 

[W]-^/2^iV^iV(0,cj2) 

when [A^] —> oo and N G Gr for some r G (0,1); this is more general than (1) 
which was obtained only for N = {n,... ,n) € Z'^. The nonfunctional version 
of LIL obtained by Wichura (1973) for any multi-parameter process with 
independent increments (in particular, for the Wiener process) allows us to 
conclude that 


limsup (2[A^] log log[A^]) = a.s. 

[A]^oo,AGGr 

We proceed now to the proof of Theorem 1.1. This is divided into several 
steps which are explained in Section 2. The remaining sections contain the 
developments that are needed to perform each step. To ease the exposition, 
we placed in the Appendix the proofs of some preliminary lemmas. 
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2. Description of the method. In this section we will indicate what are 
the main ingredients that are needed for the proof of Theorem 1.1. More 
precisely, by blending the multi-parameter blocking technique of Berkes and 
Morrow (1981) with the quantile transform technique of Csorgd and Rev&z 
(1975), we will be able to generalize to the multi-parameter case the method 
introduced by Yu (1996). 

Throughout our work we will use the letter C to denote a generic positive 
constant, independent of k. 

Let a > /9 > 1 be integers to be chosen later and no := 0. For I G Z_|_, let 


ni 


:=i:(c+/) 


2=1 


a-\-l 


-r+1 


For each k := {ki,..., kd) G Z^, we put Nk := ,..., n^^). For all k G 

Z^, we have [A^fc] ~ (a-|- l)“'^[/c]“+^. 

Let Bk := (A^fc-i,iVfc] = Us=iinks-i,nkA- Note that \Bk\ = nf=i(^? + 
k^) < We define the “big” blocks Hk and the “small” blocks Ik by 

d 

{j^ks — li ^ks — 1 T ], Ik ■ Bk\IIk- 

s=l 

Note that \Hk\ = and (2'^ — 1)[A:]^ < |7fc| < (2'^ — 1)[A:]“. We denote Uk ■= 
SiHk),Xk ■=c^‘^{Hk) and Vk := S{Ik),T^ \=a‘^{Ik). By (4), 

( 5 ) c[kr < \i < c[kr, c[kf < tI < c[kr. 


The sums over the big blocks will be used to generate a Gaussian approx¬ 
imating sequence {r]k)k which will in turn be approximated by a Wiener 
process. In order to do this, we will need an upper bound for the covariance 
of the sums over two big blocks in terms of the distance between these blocks. 
The small blocks are introduced simply to give some space between the big 
blocks, that is, to ensure that the distance between any two big blocks is 
nonzero. 

If the distribution function Fk of Uk/Xk is continuous, then one could 
use directly the quantile transform method of Csorgd and Revesz (1975) to 
approximate the variable Uk/\k by a AI(0, l)-random variable. In general, 
this assumption may not be satisfied, and, therefore one, needs to employ a 
“smoothing” technique [see Yu (1996)]. Without changing its distribution, 
we redefine the random field (wfc)fcgzd on a rich enough probability space 
together with a random field {wk)k^ 7 ,<‘- of independent random variables 

such that Wk is A^(0, r|)-distributed and {uk)k and {wk)k are independent. 
Let 


^k{uk + Wk)/, k£Z'^, 
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and Fk be the distribution function of ■ By the CLT for associated random 
fields, Fkix) <l>(x) as k ^ oo and, consequently, -Ffe(x) —> d>(x) as A: —> 
oo, where $(a;) denotes the A^(0,1) distribution function. Therefore, it is 
reasonable to consider the following A^(0, l)-random variable 

as an approximation for Let := \J+ — %)• 

In what follows we will adapt the method introduced by Berkes and Mor¬ 
row (1981) for mixing random fields to suit the special needs of an associated 
random field. 

Following page 25 of Berkes and Morrow (1981), we let r € (0,1) be ar¬ 
bitrary, p := r/8, L be the set of all indices i corresponding to the “good” 
blocks Bi C Gp, and H be the set of all points in Z'^ which fall in one of 
the good blocks. The good blocks collect those points n G which are not 
too close to the coordinate axes. From the technical point of view, their 
indices i satisfy the relationship ig > which is crucial in the proof of 

Lemma 3.9. 

To each point N € H, we associate the points AlB),..., which can be 
thought as the intersections of the hyperplanes ng = Ns,s = 1,... ,d with the 
“boundary” of the domain H; their precise definition is N^, ’ = Ng/^ \l s ^ s 
and 

:= min Ug. 

neH : n^i=N,s'^s 

Unlike the above-mentioned authors, we raise a small technical point by 
noting that H may not be a nice “L-shaped” region. This is why we consider 
the rectangles Rk := (M^, iVfc] C H, where := ((iV^^^)i,..., (iVfc'^^)d). We 
note that := {i: C C L n {i < k}. 

If U is a rectangle in and V is the rectangle in R*!, which corresponds 
to V, then we make an abuse of notation by writing W(V) instead of W(V). 



Fig. 1. 
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This convention will be used throughout this work and will occasionally 
apply to finite unions of rectangles as well. We write 


Sn = {Sn -Sn,) + S{Rk) + S{{0,Nk]\Rk), 
Wn = {Wn - Wn,) + W{Rk) + W{{Q,Nk]\Rk) 


and we use the following decomposition of S{Rk), based on the definitions 
of and Si and the fact that S{Bi) = Ui + Vi'. 


( 6 ) 


S{Rk)=J2e^+J2 





+ Vi - Y 

i^Lk i^Lk i^Lk 


In Section 3 we will show that all the sums in the above decomposition, 
except the third one, can be made sufficiently small. The third sum will 
be treated separately in Section 4 and will be approximated by W{Rk) = 

IT (Bi), via a very powerful approximation result [Theorem 5 of Berkes 
and Philipp (1979)] and a carefully chosen procedure for counting the indices 
in L. Finally, in Section 5 we will show that the terms 5((0, Nk]\Rk),W ((0, Nk]\Rk) 
can be made sufficiently small if Nk € Gt, and the differences — — 

ITvfe are small if G Gp. This will conclude the proof of Theorem 1.1. 


3. The “good” blocks. In this section we will show that all the sums in 
the decomposition (6) of Sr^ , except the third one, can be made sufficiently 
small. 

In order to treat the first sum of this decomposition, we need to evaluate 
the precision of the approximation of ^k by Vk- This will be given by the 
speed of convergence in the CLT. In this paper we decided to use the result 
obtained by Bulinski (1995), under the assumption that the covariance co¬ 
efficient u{n) decays exponentially as n ^ oo. Under this assumption, this is 
the sharpest speed of convergence in the CLT when d = 1, s = 3 [see Birkel 
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(1988b)]. We note in passing that in the case d = 1, a different speed of 
convergence in the CLT was developed and used by Yu (1996) for associated 
sequences with a power decay rate of the covariance coefficient; however, 
the exponential decay rate of u{n) was eventually needed for the strong 
invariance principle. The problem of whether or not the strong invariance 
principle continues to hold for associated random fields with a power decay 
rate of covariances is still open even in the case d=l, and we do not attempt 
to tackle it here. 

Lemma 3.1 [Theorems 1 and 2 of Bulinski (1995)]. Suppose that (Cl) and (C2) 
hold and let s := 2 + r + 6. Then for any finite subset V C Z'^, 

sup \Fv{x) - <l>(a;)| 

^ f ClYj • (ct 2(Y))-*/2 • (log(|Y| + tfs<3, 

- \ ClYj • (a2(Y))-3/2 . (iog(|Y| + i))-^, ^fs>3. 

The next result is a generalization of Lemma 3.2 of Yu (1996) to the case 
d>2, in the case of an exponential decay rate of u{n). Its proof is routine 
and is given in the Appendix. 

Lemma 3.2. If (Cl) and (C2) hold and 2rQrj{2 + r) < affi < 2(1 + 
r )/(2 + r) with ro := max{l, (r + 5)“^}, then 

sup |Tfc( 3 :) - ^>(x)| < and sup \fk{x) - f{x)\ <C, 

where fk{x) is the density function of f,k fix) is the Y( 0 , 1 ) density 
function. 

Using Lemma 3.2 and an argument that was introduced in the proof of 
Lemma 3 of Csorgo and Rev&z (1975), we get the precision of the approxi¬ 
mation of by %. 

Lemma 3.3. Under (Cl) and (C2), we have 

\^-\Fkix)) -X\< C[k]-^^d/(2+r)-Ky2} ^ 

provided that jxj < K^J\og[k], where 0 < 77 < ^2rj3ji2 + r). 

Next we give the precision of the approximation of by r]k in terms of 
the L^-distance. For this we will need the following lemma which gives an 
upper bound for the moments of order 2 -|- r, generalizing an older result 
of Birkel (1988c) in the case d = 1. In particular, this lemma shows that 
iXj)j^ 2 f finite r-susceptibility (as defined in the Introduction). 




INVARIANCE PRINCIPLE FOR RANDOM FIELDS 


9 


Lemma 3.4 [Corollary 1 of Bulinski (1993)]. Suppose that (Cl) and 
(C2') hold with v > dvo, where vq := r(2 + r + 5)/{25) < {d — 2)~^ if d>3. 
Then for any V € A, 

E\S{V)\‘^+^ <C\V\^+^/^. 

Using (5), Lemmas 3.3 and 3.4, and employing the same technique that 
was used in the proof of Lemma 3.10 of Yu (1996), we get the following 
result. 

Lemma 3.5. Under (Cl) and (C2), we have 
where Eq := 2r^/3/{(2 + r)(4 + 3r)}. 

The next result will show us that the first sum in the decomposition ( 6 ) 
of S{Rk) is small. 

Lemma 3.6. Suppose that (Cl) and (C2) hold and (5 > (l + 2/r)(3 + 
4/r). Then there exists ei > 0 such that for every k G with 0, 

\ei\<C[Nk]^/^~^^ a.s. 

i&Lf. 


Proof. Let <7 > 0 be such that a — Eq + 1 <2q < a — 1 (this is possible 
since Eo> 2 hy our choice of /3). By Chebyshev’s inequality and Lemma 3.5, 
we have 

P(|ei| > [f]^) < [i]-{29-(“-eo)} Vz G Z^. 

By the Borel-Cantelli lemma, it follows that \ei\ < C'[i]'?, Vi G Z'^ a.s. and, 
hence, le^l < < C[k]'i^^ < < C[iVfc]V 2 -^i a.s., 

where 0 < e[ < {a — 1)/2 — q and Ei := E'il{a + 1). □ 

The proof of the following lemma is given in the Appendix. 

Lemma 3.7. If {C2') holds with d<v < 2d, then 

(7) = 

where V is a finite union of rectangles in A and 5 q'.= v jd—Y. 

Remark. Relationship (7) is exactly Dabrowski’s (1985) condition for 
the FLIL for associated sequences. 
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Lemma 3.8. Suppose that (C2') hold with d < v <2d and (3 > 3/(5o, 
where (5o '■=v jd — 1. Then for every k € Z'^ with 7 ^ 0, 


ieLj. V V I *1 / 


a.s., 


where ao := l/{ 2 (a + 1 )}. 


Proof. Note that Ui := a — (A? + Tf)/\Bi\ > 0, by (4) and the associ¬ 

ation property. Using (7), we have 


\^i 


a^- 


A? 


\Bi\ V m 


+ 


1^*1 

\Bi 


a 



and, hence, by Chebyshev’s inequality. 

By the Borel-Cantelli lemma, it follows that y^\Bi\air]i < Vi G Z^j[_ 

a.s. and, hence, T,i&Lk V\W\aiVi < 

a.s. since [A:] ~ (a -|- . □ 


The hnal result of this section shows that the last two sums in the de¬ 
composition (6) of S{Rk) are small. 

Lemma 3.9. If a — P > 2 + 4/ p, then for every k G Z'^ with Lk 7 ^ 0, we 
have 

Y,h\<C[Nk]^/^-^° a.s. and 1^*1 ^ 

Proof. For the first inequality, we follow the proof of Lemma 8 of 
Berkes and Morrow (1981). Note that T = Uf=i where Ii{s) are disjoint 
rectangles with |/i(s)| < Ci^Yis'^s^'i'- Hence, Vi = with vps) := 

J2jGii{s) ^j- 

By Chebyshev’s inequality and (4), 

p(b*(s)i > < c[i]-(“-2)|i,(s)i < n e 

s'^s 

< •-(a-/3-2-2/p) < ^j.p(a-/3-2-2/p)p/2 

for every i £ L^. (As in the proof of the above-mentioned lemma, we used 
the fact that i £ implies that ig > CHs' and, consequently, ig > 
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Since {a — (3 — 2 — 2 // 3 )p /2 > 1 , the result follows by the Borel- 
Cantelli lemma. A similar argument applies to Wi, since E{wf) = rf < 

C\h\ = CY.i=i\h{s)\. 

□ 

4. The approximation theorem. In this section we will verify that the 
third sum in the decomposition ( 6 ) of S{Rk) can be approximated by W(Rk), 
where W is a d-parameter Wiener process with variance a^. Some prelimi¬ 
nary lemmas are needed. 

The next result follows exactly as Theorem 2.1 of Yu (1996), using Lemma 3.2. 

Lemma 4.1. If (Cl) and (C2) hold and 2rQr jifl + r) < a/P < 2(1 + 
r)/( 2 -|-r) with tq := maK{l,{r + then for any 0 < 6 < 1/2 and all 

Eimh,) < c{([i][j])-“/2L;(u,u,)}®/(i+^). 

The next lemma gives a generalization of relationship (3.11) of Yu (1996) 
to the multi-parameter case. 

Lemma 4.2. If (C2) holds, then 

E{uiUj) < 

where Mij := max^, j) - 1 ) and Ms{i,j) := max(i 5 , j^), s = 1 ,..., d. 

Proof. Let d := mink^Hi d{k,Hj) be the distance between Hi and Hj, 
where d{k, Hj) := mink'^Hj \\k — k'\\. Then dk ■= d{k, Hj) — d>0 y k £ Hi, 

E{uiU,) = E E E{XkXkP < E u{d + dk) < Ce-^^ E ^ 

k£Hik'£Hj keHi k&Hi 

andd = m.ayis=i,...4mm.k£Hi,k'£Hj \ks-k'P = max^, -L+ 

id)} > where nig = ms{i,j) := mm{is,js) and Mg = Mg{i,j). □ 

In order to prove our approximation theorem, we need to be able to 
“count” properly the indices in L, that is, to dehne a bijection fj : Z 4 . ^ L 
satisfying certain properties. This will be given by the following lemma, 
whose proof can be found in the Appendix. 

Lemma 4.3. There exists a bijection if : Z+ ^ L such that 

(8) l<m 3s* = s*{l,m) such that'ijj{l)g* <‘ip{m)g* 

(9) 3mQ £ such that m < C[p{m)p° Vm>mo 
for any 70 > (1 -L l/p){l - 1 /d). 
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We are now able to prove the desired approximation theorem. 


Theorem 4.4. Suppose that (Cl) and (C2) hold, a > 3(1 + 1/ p){l — 
1 /d); (3 > ( 2 /p)(l + l/p){l - 1 /d) and 2 ror /(2 + r) < a/fd < 2(1 + r )/(2 + r) 
with ro := max{l, (r + d)“^}. Then without changing its distribution, we 
can redefine the random field {Xj)j^zf ® enough probability space 

together with a d-parameter Wiener process W = {Wpt G [ 0 , 00 )*^) with vari¬ 
ance , such that for every k G Z/. with Lkfi^0, 




Vi- 


W{Bi) 

cr\/\B~\ 




a.s., 


where ao := 1 /{ 2(1 + a)}- 


Proof. Let 0 < d < 1/2 be such that a{{\ + l/p)(l — l/d)}“^ > 1 + 1/d 
and choose 70 such that (1 + l/p)(l — 1/d) < 70 < min{ad/(l -\- d),j3p/2}. 
Let : Z_|_ —> L be the bijection given by Lemma 4.3. 

We will apply Theorem 5 of Berkes and Philipp (1979) to the sequence 
Ym := £ Z-|_ of random variables and the probability distributions 

Gm ■= 1V(0,1), m G Z-|_ and we will prove that for each m G Z+, m > 2 there 
exists some pm > 0 such that 


( 10 ) 


E exp 



E exp 


m —1 

i E 


i=i 


EeKpfiitmYm] 


El Pm 


for all ti,... ,tm with E U^, where Um > lO^m^. 

Then, by the above-mentioned theorem, without changing its distribu¬ 
tion, we can redefine the sequence (Tm)mez+ on a rich enough probability 
space together with a sequence {Zm)mez+ of independent A^(0, l)-random 
variables such that 


-P(|hm Z,ji\ ^ CTm) — (^m V?7i G Z_|_, 

where am E C{l7m^''^logUm + exp(-3Um^/W)mE^Um^ + pYlEm'^^''^}- We 
will prove next that 

( 11 ) am E Cm~‘^ for m large. 


Then, by the Borel-Cantelli lemma, \Ym — Zm\ E Cam, Z+ a.s. Us¬ 

ing a straightforward d-parameter generalization of Lemma 4 of Csorgo and 
Revesz (1975), without changing its distribution, we can redefine the se¬ 
quence {Zm)mez+ on a richer probability space together with a d-parameter 

Wiener process with variance such that Zm = W{B^(^m)) / {^ I ), Vm G 
Z+. Hence, 


Vi 


W{Bi) 




E Ca^ 


Vi £ L a.s. 
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and because \Bi\ < \Bk\ < CfA:]", \/i^ and I]iGZ+ 


i^L]z 


m- 


W{Bi) 

(y\/\Bl\ 


< c[kr'^ Y. «b-Hd ^ ^ c[iVfc]i/2-«o 


We proceed next to the verification of (10) and (11). By Lemmas 4.1 and 4.2, 
we have 

E{YiYrn) < C{([V^(/)][^(m)])-“/2E(u,^(;)U^(^))f 


[For the last inequality above we used (8) to obtain an s* = s*{l,m), for 
which Ms*{'ilj{l),'tp{m)) ='ip{m)s*', since ip{m) G L, we have > 

'ip{m)s* — 1 > C['ip{m)Y^‘^.] By Lemma 2.2 of Dabrowski and Dehling (1988), 
the left-hand side of (10) is smaller than 2X)[(L7^ \titm\E{YiYm), which is, in 
turn, smaller than 

m—1 

{ m—l 

1=1 


+ (m 


l)tl^[ip{m)] 


-ae/{i+e) 


< ^^-xemmrf>/Vii+e)j2tf < Ce-^emm)]^^/Vii+9)ul ~ 

1=1 

for m large enough. (In the second inequality above, we used the fact that 
m < C['ilj{m)]'^^E^+9) ^ which follows from Lemma 4.3 by our choice of jq.) 
Finally, relationship (11) follows if we take Bm '■= with q> 8. Clearly, 
log Um < rn~‘^ and exp(—3[/m ^/16)m^/^[/m < exp(—2[/rn ^/16) < 
for m large enough. We have 

^V2f/m+l/4 ^ ^-Xemm)fP/y{2+2e)^q{m+5/A) < ^-2 

since {2-|-(?(m-|-5/4)} logm < < C'[V’(m)](^+^)'>'° < C['i/i{m)]^9/‘^^ for 

m large enough. This concludes the proof of the theorem. □ 


Remark. A similar argument can be used to give a simplified proof 
for Theorem 2.5 of Yu (1996) (in the case d = 1). More precisely, one can 







14 


R. M. BALAN 


check directly the condition of Theorem 5 of Berkes and Philipp (1979) for 
the sequence (%)a:>i of random variables and the probability distributions 
Gk = fV(0,1), /c > 1 (as we did above). We obtain in this manner a sequence 
{Zk)k>i of independent A^(0, l)-random variables with P(|r/fc — Zk\> ak) < 
ak and ak < Ck~‘^. Without changing its distribution, we can redefine the 
sequence {Zk)k>i on a richer probability space together with a standard 

Brownian motion W = {Wt',t G [0, oo)} such that Zk = W{Hk)/\J+ 

where Hk := (14-i, 14] and 14 := Z]i=i(^i + T'i)- Since Af + rf < 

for i < A; and X]i>i Oi < oo, this gives immediately the desired approximation 


k 

E 

i=l 




2 + t 2 

I ' ' I 


Vi 


W{Hi) 

\/E+E 


k 

i=l 


a.s. 


5. The remaining terms. In this section we show that the termsS'((0, Nk]\Rk), W ((0, Nk]\Rk), Sn 
S'tVj,, Wn — Wn^ can be made sufficiently small if G G,-. 

Note that {0,Nk]\Rk = If we let Ds{N) :=max„<^(s) \Sn\ 

and Ds{N) := |114,|) for each s = 1,... ,d and N & H, then 


d 

Si{0,Nk]\Rk) 

S = 1 


d 

Wi{0,Nk]\Rk) <Y.‘^'^~"Ds{Nk). 

S = 1 


On the other hand, (0, A^]\(0, N'fe] = wherex 

nsej=(0) ^ks] and the union is taken over all nonempty subsets J of {1,, d}. 
Let := max |S'(/^'^^)| and := sup |1T(/^'^^)|, where the maximum 

and the supremum are taken over all N with Uk^ < Ng < Uk^+i, V s & J. We 
have 


max \Sj\f — Sn^ I < , 

Nk<N<N^+^' k ^ 


sup |lP7V-fLwJ< 

Nk<N<Nk+i j 


We note in passing that the arguments that are valid for the terms 
depending on the original random field {Xj)j^z'^ f*® applied to the 

terms depending on the Wiener process IT, since 1T(T) =J2j£V^j^ 4T G 
A, where Xj := W{{j — l,j]) are independent ^"(0, cr^)-random variables. 
Clearly, {Xj)j^ 2 ,d is a weakly stationary associated random field with zero 
mean and covariance coefficient u{n) = 0, Vn > 1. 


Lemma 5.1. (a) Suppose that (Cl) and {C2') hold with v > dvQ and 

z^o := r{2 + r + 6)/{25) < {d — 2)“^ if d> 3. Then there exists xq such that 
VT G .4, Vx > xo, 


P{M{V) > x\V\^G) < Cx-^^+O 
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where M{V) := max{|S'(Q)|; Q C y,Q £ -^}- 

(b) If (Cl) and (C 2 ) hold, then there exists 7 > 0 such that ^ A, 

p{M{v) > |y|^/^(iog |y|)'^+^) < c|yr^, 

where M{{a, b]) := max{|5(Q)|; Q = {a,c],a < c <b}. 

Proof, (a) Using Lemma 1 of Bulinski and Keane (1996), the Markov 
inequality and Lemma 3.4, we have P(M(U) >a:|U|^/^) <2P(|5(U)| >x\V\^/‘^/2) 

Cx-( 2 +U|y|-(l+U 2 )^| 5 (y)| 2 +r < fj^-{2+r) ^ 

(b) This follows exactly as the second inequality of Lemma 7 of Berkes and 
Morrow (1981), using the moment inequality given by Lemma 3.4 and the 
rate of convergence in the CLT given by Lemma 3.1. This rate is sharper than 
the rate of Lemma 5 of Berkes and Morrow (1981). To see this, we use (4) and 
we note that sup^-gj^ \Fv{x) — 4>(x)| is either smaller than 
if s < 3, or smaller than if s > 3; in both cases a suitable choice 

of e > 0 gives us the rate ClUl"* for some t G (0,1). We also note that the 
requirement |U| G G,- is not needed. □ 

The next result follows exactly as Lemma 6 of Berkes and Morrow (1981), 
using Lemma 5.1(a). 

Lemma 5.2. If a> 16/(3r) — 1, then 

max Ds{Nk) < a.s., 

s=l,...,d 

max bs{Nk)<C[Nk]^/‘^-^ a.s. 

s=l,...,d 

for every Nk G and 0 <e < r/32. 

The following result follows exactly as Lemma 9 of Berkes and Morrow 
(1981), using Lemma 5.1(b). 

Lemma 5.3. Let 7 be the constant given by Lemma 5.1(b). If ot> 2 / 7 , 
then 

maxM^'^^ < a.s., 

maxM^'^^ < C[Nk]^/^-^ 
for every G Gp and 0 <e < p/ ( 80 :). 


a.s. 
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APPENDIX 


Proof of Lemma 3.2. Using Lemma 3.1 for V = Hk and relationship 
(5), we obtain that |Ffc(a^) —‘h(x)| is either smaller than 

if s < 3, or smaller than if s > 3. If a/(3 > 2rQrl{2 + r), then a 

suitable choice of e > 0 allows us to conclude that |Ffc(x) —‘h(3:)| < Vx G 

R. The first inequality follows by a change of variables. 

For the second inequality we use a technique similar to that used to 
prove relationship (3.3) of Yu (1996). Let ipk{t) ■= E[exp{it^k)],‘fk{'t) ■= 
E[exp{ituk/\k)] and (p{t) =exp(—1^/2). Since (A| + t|)/A| < C, we have, 
for any T > 0, 

\fk{x) - fix)\<-^ f \(pk{t)-(p{t)\dt 
^'TT J— oo 


c r°° ( 

< 7 iz/ \^kit)-(p{t)\exp}^-^ 


< 


27T j —oo 

c 


ds 


2tt 


2T[k] 


-r/3/(2+r) 


+ 


c 


vr J|t|>T 


exp 


2Ai 


< C . r[tl-''’/P+’'> + ^ . ^exp|-|T|. 


dt 


Since A|/r| < the conclusion follows by choosing T = C[k]'^ with 

a — P < q < rP/{2 + r). Such a choice is possible if a//3 < 2(1 + r)/(2 + r). 

□ 


Proof of Lemma 3.7. First we claim that it is enough to prove (7) 
for “squares,” that is, for rectangles V = (m,re] € A for which — rrig = 
I, V s = 1,..., d. To see this, we note that each rectangle V can be written as 
a hnite union of disjoint squares: V = UlLi By the association property 
o'^(^) > Ef=i<7^(Ui) and 




1 


1^1 


1^1 S'"*' 


a2- 


^Hv^) 

1^*1 




because 0 < do < 1- Let us now prove relationship (7) for a square V = (m, n] 
with Us — nis = 1, Vs = 1,..., d. Note that \V\ = 1^. By stationarity. 


cj2(U) = |U| .r(0) + ^ 

— (n—m—l)<2<n—m—1,27^0 5=1 

= |U|- E (-1)'^'"^ E c{Kp)-r{i), 

||2||</ —1 07^i^ C{l,...,ci} ||2||</ —1,27^0 
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where c{K, i) := K* I • Since cr^ - E||i||<z-i ^(0 = E||i||>z ^(0 = u{l) 

and c{K,i) < \ V\ if ||i|| < / — 1, we have 

Y. M ^ ciK,i)-rii) 

0^A'C{l,.,.,d},lA'|odd ' ' 

<c\vr/'^+ Y E 

0^KC{l^...,d},\K\ odd || 2 II</— 1 , 27^0 

We used the fact that u{l) < Cl~'^ = C\V\~'^/'^ and r{i) < u(||i||) < < 

for any i e Z'^, where [i] = l^sl- 

Proof of Lemma 4.3. The idea of the proof is based on the following 
simple observation in the case d = 2. For each m G Z_|_, m > 2 with (m, m) € 

L, there exists a /c|(m) > m such that (/ci,m), (m, fci) G L for every m < 
ki < kKm). Therefore, to each vertex {m,m) G L, one can associate an “L- 
shaped” region L{m) consisting of 2{kl(m) — m} + 1 points in L. In view of 
the desired property (8), we will count consecutively the indices in L(2), L(3), 
and so on. To verify property (9), we note that k G L{m) implies [k] > m?. 

We begin now the proof for arbitrary d>2. Let m G Z_|_,m > 2 be such 
that (m,..., m) G L and k = {ki,..., k^-i, m) G L be such that kg > m, Vs < 
d. This implies that all the vertices of Bj. are in Gp, and, in particular, rim > 
y s < d. Since m is fixed, this cannot happen for infinitely many kg’s. It 
follows that for each s = l,...,d—1 , there exists a k*{m) > m such that kg < 
kg{m). We note that k*{m) < if m is large enough. This argument 

shows us that we have a maximum number of k*{m) : = I\i=l{k*g{'m) - m} 
points of the form {ki,..., kd-i,m) in L, with kg > m, Vs < d. 

By symmetry, we can repeat this argument for each of the axes. We let 
Lg{m) := {A: = (fci,..., /c^-i, m, kg,..., kd-i)',m < kg/ < k*,{m), Vs' < d} for 
every s = 1,... ,d. The “L-shaped” region corresponding to the index m is 

d 

L{m) := y Lg{m) U {(m,... ,m)}. 

S = 1 

Note that |L(m)| = dk*{m) + 1 and that k G L{m) implies [/c] > rrd. Clearly, 

L C IJ^ L{m) [note that in the case d = 2, we actually have L = IJ^ L(m)]. 

Next we count consecutively the indices in L(2),L(3), and so on, that is, 
we define a bijection ip: Z+ ^ IJ^ L{m) such that Vz G Z+, 

m—1 ra 

Y \L{i)\<z<Y\m\ 

1=2 1=2 


ip{z) G L{m). 
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The bijection ip clearly satisfies condition ( 8 ). To verify (9), we note that 

m d—1 d—1 m 

z<d'^Yl{k*{l)-l) + m<dY[ -l) + 'm 


1=2 s=l 
d—l / m 


s=l 1=2 
\ d-1 


d—1 / m 




s=l \ 1=2 


. 1=2 


s=l \ 1=2 


< (^^(i+VpKrf-i) < < C[<p>{z)\^° 


for m large enough and 70 > (1 + 1/p) (1 — l/d) arbitrary. Finally, define the 
bijection 1 /: Z+ ^ L such that < (p~^{k), Vk £ L. The result follows 

since if zi,Z 2 £ Z+ are such that V’(- 2 i) = +( 2 ^ 2 )) then zi < Z 2 - □ 
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